Full Counting Statistics of Cooper Pair Shuttling 
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The Cooper pair shuttle is a simple model system that combines features of coherent and incoher- 
ent transport. We evaluate the full counting statistics (PCS) of charge transfer via the shuttle in the 
incoherent regime. We describe two limiting cases when the PCS allows for classical interpretation. 
Generally, the classical interpretation fails yielding negative and imaginary "probabilities". This 
signals that superconducting coherence survives even in incoherent regime. We evaluate the current 
noise in some detail. 

PACS numbers: 



The Josephson cffecti consists in coherent Cooper pair 
transfer through a tunnel junction between two bulk 
superconductors. It results in a dissipationless current 
present as a ground state property of the system. This is 
in contrast to dissipative electron transfer in normal con- 
ductors. In latter case, it is possible to evaluate the full 
counting statistics (FCS) of charge transfers— in terms of 
classical probabilities. Thereby, the transport properties 
are fully understood and characterized: one can predict 
not only the average current, but also the current noise 
and all higher moments of the current distribution func- 
tion^. 

One can also access the FCS of a superconducting 
Josephson junction^. Due to gauge symmetry breaking 
in the superconducting state, this FCS can not be in- 
terpreted in classical terms yielding negative "probabili- 
ties" . The use and the interpretation of the FCS in this 
case is that it determines the quantum evolution of the 
detector that measures the currenii^. 

A recent proposal puts forward an interesting way to 
transfer Cooper pairs: to shuttle them (controUably) be- 
tween the superconducting electrodes^. Although the 
original proposal puts emphasis on mechanical degrees of 
freedom, a Cooper pair shuttle is essentially a supercon- 
ducting single electron transistor (SSET) with variable 
time-dependent Josephson coupling to the superconduct- 
ing leads^. The Coulomb island of this SSET is brought 
in contact with one electrode at a time, the electrode 
contacted being periodically altered^. 

It is interesting that the supercurrent between the 
leads is achieved as a result of a non-equilibrium driven 
process, generally accompanied by dissipation. The co- 
herence of Cooper pairs transferred is determined by the 
coherence of different charge states in the island, this be- 
ing mostly affected by fluctuations of the gate voltage^. 
The transport results from the interplay between the co- 
herent Josephson coupling and decoherence effects. This 
is why the Cooper pair shuttle presents a model to bridge 
between the limits of coherent and incoherent transport. 
The model is simple indeed, one can restrict the consid- 
eration to just two quantum states. 



Thus motivated, we have analyzed the FCS of the 
Cooper pair shuttle in the most interesting incoherent 
regime, where no net supercurrent is shuttled. The 
regime is achieved in the limit when the voltage fluctua- 
tions are classical. If the fluctuations come from an en- 
vironment, this implies that the temperature of the envi- 
ronment exceeds the relevant energy scales of the shuttle. 
Naively, from the absence of the net current one would 
conclude that no charge transfer occurs in the system. 
However, transfers do occur, the current is zero only in 
average, and the FCS presents a convenient way to reveal 
this circumstance. This stipulates the understanding of 
transport properties of the shuttle and, generally, the in- 
terplay between the coherent and incoherent transport. 

The results are as follows. If the period of the shuttling 
is sufficiently long for decoherence to be accomplished, 
the FCS can be interpreted in terms of classical elemen- 
tary events: Cooper pair transfers. During the shuttling 
cycle, either no transfer takes place or a pair is trans- 
ferred in either direction. There is an apparent similarity 
with the FCS of the pumping in normal systems studied 
in^J^. The FCS in the opposite limit of short cycles al- 
lows for alternative classical interpretation in terms of 
a superconducting current randomly switching between 
two opposite values. 

In the general intermediate situation, the FCS can not 
be interpreted in classical terms. An attempt to evalu- 
ate the probabilities per cycle yields negative and even 
imaginary values. This is a clear signature of the fact that 
the superconducting coherence survives strong dephasing 
although this coherence does not manifest itself in net su- 
perconducting current. It was recently explainedii that 
any FCS can be characterized directly so one would not 
have to measure higher cumulants of the current noise 
one by one. However, the immediate physical value mea- 
sured would also depend on the properties of the concrete 
detector. This is also a way to experimentally observe the 
FCS of the Cooper pair shuttle. 

The system of interest is presented schematically in 
Fig. [TJ It consists of a Cooper pair box, or Coulomb 
island, connected with Josephson junctions to the super- 
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conducting leads L and R. Conventionallyf«iiiiS we as- 
sume the separation of energy scales, A » Ec » Ej, 
A being the superconducting energy gap, Ec being the 
charging energy of the box, and E,j being a typical 
Josephson junction energy. Under these conditions quasi- 
particles degrees of freedom are not involved in the sys- 
tem's dynamics. In addition, the gate voltage is chosen 
to bring two charge states of the box, say, |0 > and |1 > 
to the degeneracy point. Under these assumption, we can 
restrict the consideration to these two states only and the 
shuttle is described by the Hamiltonian 



H{t) = 2eV{t)a^ 
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which we write in terms of 2 x 2 Pauli matrices a^, 
a± = {axi<Ty)/2 in the space spanned by |0 >, |1 >. V{t) 
is proportional to the deviation of the gate voltage from 
the value corresponding to the degeneracy point and (f>R, 
(pL are the phases of the macroscopic superconductors. 
Only (j) = (fiR — (f>L is a physical quantity which we will 
assume to be fixed; this can be obtained by closing the 
circuit endpoints in a loop pierced by a constant magnetic 
flux. Ej^'^\t) are time-dependent Josephson energies of 
left and right junction. For the sake of concreteness, we 
assume stepwise periodic variation of both Josephson en- 
ergies with time as shown in the lower panel of Fig. 1. 
Each lead is contacted during time interval tj and the 
box contacts no lead during time interval tc for each pe- 
riod T = 2(tj + tc)- We call tj and tc the "Josephson 
contact" and the "free evolution" time respectively. An 
idea to realize the required time dependence in Joseph- 
son coupling by means of an experimentally realizable 
SQUIDs based device has been proposed in Ref. 0. 

We assume V{t) to be a classical stochastic vari- 
able with white noise statistics, (T^(i)}stoc ~ ^ ^'^'^ 
(^(i)^(i'))stoc. = lhy{2emt-t'), where defines 
the average over the fluctuations. Thus defined, 7 is just 
the inverse decoherence time of the two charge states. If 
we neglected the fluctuations, the time evolution of the 
system would be fully coherent governed by the coherent 

part of the Hamiltonian He = \ H) .In this case the 

\ / stoc. 

quantum state of the central grain would acquire dynam- 
ical phases 29 = Ejtj/h during the Josephson contact 
times and 2x = 2eVotc/fi during the free evolution times 
(we assume that Vq is only present during both intervals 
of the free evolution time) . 

With fluctuations, the shuttle will be described by a 2 x 
2 density matrix that obeys the following Bloch equation: 



H,it)p~ pH,(t)] ~2j{p~a,pa,) . (2) 



dp^ 
dt 

The only stationary solution of this equation is quite 
trivial: p oc 1, this corresponds to the absence of the av- 
eraged superconducting current. Therefore, the shuttle 
operates in incoherent regime. This is a combined effect 
of the decoherence term and Josephson coupling. In the 
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Figure 1: Upper panel. Cooper pair shuttle consists of a 
Cooper pair box coupled to superconducting leads through 
Josephson junctions. Two charge states are tuned to degen- 
eracy point by the gate voltage. Lower panel. The specific 



time dependence of Ej 
the shuttling. 
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within a single period T provides 



absence of Josephson coupling, voltage fluctuations can 
not cause transitions between the charge states so no re- 
laxation takes place. With Josephson coupling, the volt- 
age fluctuations cause transitions between the stationary 
states separated by energy 2Ej at Vq = 0. Classical 
voltage fluctuations result in equal transition rates with 
increasing and decreasing energy. In fact, the ratio of 
these rates is given by Bolzman factor cxp{2Ej /Tf,), T\j 
being the temperature of the environment producing the 
fluctuations. As shown in Ref. 0, proper account of this 
factor leads to anisotropic p ^ \ and to the net supercur- 
rent vanishing at Tf, ^ Ej. This defines the incoherent 
regime under consideration. 

Let us evaluate the FCS. One starts with the assump- 
tion that a transport process can be characterized by the 
probabilities Pt{N) of N electrons transferred through 
the contact in a time interval t and attempts to compute 
the characteristic function of this probability distribution 
defined as^ 



=.-e(A,r) 



XN 



(3) 
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The general quantum expression for this function can 
be obtained by coupling the system to a detector for a 
measuring time r and interpreting the detector readout 
in terms of charge transfer—. The FCS in this case is 
defined as integral kernel that relates the initial and final 
density matrices of the detector, and reads 



Tr 



W+A(T,0)/5(0)iYl,(T,0) 



(4) 



p'(r) 
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where Wa(t, 0) = T Hx{s)ds jg ^-^^ unitary evolu- 
tion operator corresponding to the modified Hamiltonian 
Hx = H + XI /e, I being the operator of the electric cur- 
rent. Two such operators provide non-unitary evolution 
of the initial density matrix of the system p{0) into p' {t). 

For our shuttle model, the implementation of this gen- 
eral scheme is especially simple since the density matrix 
is just 2x2 matrix. Following Ref. 0, we also gauge 
the counting field A to the left electrode, this yields 
H\ = Hc{<t>L 4>L + A). We derive an equation for 
p' that looks very similar to the Bloch equation ^ 

% = [h+xp' p'H^x) - 27 [p' - a.p'a.) . (5) 

The difference is that the Hamiltonians governing evolu- 
tion of "bra"'s and "kef's differ by opposite A shifts. 

This system of 4 linear equations can be solved with 
initial conditions p'{to) to obtain a linear map M that 
gives p' after a cycle, p'{to+T) ~ Alp' {to) All four eigen- 
values /ii of M satisfy the condition \p,i\ e [0, 1]. Since 
we study statistics of low-frequency fluctuations, t T 
by definition. Therefore, the FCS is determined by the 
eigenvalue with the greatest magnitude. 



-6(A,r) 



^{Xy^"^ , where 1— 
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t=l,...A 



This is in accordance with the method and the resulti^ for 
FCS of charge transport described by a master equation. 

The most transparent way to present the FCS is to 
define the probabilities pn to transfer N electrons per 
cycle, as it has been done in Ref. to characterize the 
pumping of normal electrons. 



Pn 
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/2(A). 
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The FCS in the limit considered is always an even func- 
tion of A, so that the pairs are shuttled in either direction 
with equal probability, pn ~ P^n- 

We analyze the FCS in several limiting cases. First 
we consider the limit of long periods l/tcA/ij ^ 7 ^ 
Ej/h. This is in fact an adiabatic limit since the relax- 
ation and decohcrencc arc fully developed within each 
time interval tc,tj. In this case, 



li{X) — cos^ A 
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so that, each shuttling between the superconductors 
transfers either one Cooper pair or none, this occurs 
with equal probabilities. The pair is transferred with 
equal probabilities in either direction. It is interesting to 
note that this simple result is quite general and relays 
on neither the periodicity of shuttling nor the concrete 
time dependence of Ej{t) provided the adiabaticity is 
preserved. The shuttle may even return (several times) 
to the same superconducting terminal before contacting 
the opposite one: the charge transfer is in this case as- 
sociated with two trips between the opposite terminals. 



Leading corrections to adiabatic FCS are exponentially 
small, ~ exp(— 2tj7). 

Another interesting limit is that of small Joseph- 
son couplings Ej ^ and sufhciently long cycles 
h/tc,fi/tj ^ Ej. In this case, the charge relaxation 
time ~ jfi?/Ej may be long, exceeding both decoher- 
ence time and cycle duration. The relevant parameter 
exTp{—tjEj/{h'^j)) = f determines the efficiency of the 
charge relaxation during a cycle. The FCS becomes more 
complicated, for / ~ 1 it is given by 

iiW = kl + f) ms' A + /sill" \ + (l-f)x 



cos^ A 1 



1 + / 



^ sin^ A 
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and, generally speaking, all pn ^ 0. This is because the 
charge can not completely relax during a cycle providing 
"memory effect" so that charge transfers in different cy- 
cles are not independent and elementary event of charge 
transfer can encompass several cycles. Still, probabili- 
ties remain positively defined owing to the fact that in 
the limit of long cycles the FCS depends neither on the 
superconducting phase cj) nor on dynamical phases 9, x- 
Beyond this limit, the FCS does depend on <f). As men- 
tioned in Ref. 0, the classical interpretation in this case 
may fail and one can not assure that p„ are positive or 
even real. Indeed, we have found this in numerical cal- 
culations at jT c± 1 (Fig. This clearly signals that 
the superconducting coherence remain in the system al- 
though the decoherence completely destroys the average 
supercurrent. 

This becomes evident when analyzing the opposite lim- 
iting case of very short cycles ■jT <C 1. In this case, 
the cumulants of the current arc contributed by FCS at 
A ~ ^ 1. In this region, the FCS reads 



T 




(10) 



where Is is the superconducting current in the absence 
of decoherence 7 = and 7 is the relaxation time in the 
limit 7^0: 



2e da 



tr 



-Pr 



where a = arccos (cos^ 6* cos 2% — sin^ 6* cos 0) and Pj, 
Pc are always positive polynomial functions in sin(-) and 
cos(-) of the phases 0, x, Q- 

Pci4>,x,0) = 2sin2(26i)(l + cos(/)COs(2x)) 
Pj{(j),X,6) = 2 (1 + sin2 6lcos(/icos(2x)) + 
- cos(26l) (cos2 6icos2(2x) 



sin^ 6 cos^ I 



FCS of this type corresponds^ to a supercurrent ran- 
domly switching at time scale I/7 ^ T between the op- 
posite values ±/s, so that the coherence is preserved at 
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this time scale. The charges are transferred in long trains 
of ~ 1/7^ of elementary charge, this leads to significant 
low- frequency current noise. 

Low-frequency noise S{uj = 0) is obtained di- 
rectly from the definition of the FCS, 5'(0) = 
-e2/riimA^o5^(M/5A2). From the Eq. (O we obtain 
5(0) = 2eVrx in the adiabatic limit. The 

noise is enhanced in the opposite limit of short cycles, 
S'(O) = 2/g/7 » e^/T and is sensitive to all dynamical 
phases. The numerical results in the intermediate regime 
(Fig. [3]) show in addition a quasi-oscillatory dependence 
on the dynamical phase 6. 




Figure 2: Fourier components, pjv, with = 0, 2, 4 of /i(A) 
plotted versus Ej /fi'y. Each line is labeled by the correspond- 
ing value of A^. Adiabatic regime (a): ')tj = 5; positive 
probabilities. General situation (b-c): 'ytj — 0.9; pjv can be 
negative and even immaginary. Other parameters are fixed 
to tj = tc, X — ■'■"/S) 4> = ■'i'/4 for all plots. 



Figure 3: Low- frequency noise in the intermediate regime for 
difi'erent values of 'ytj, as expressed by each line label. The 
noise and its oscillations with increase with decreasing "/tj. 
Other parameters are fixed to tj = tc, X = ^/5, = 7r/3 for 
all plots. 
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